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ABSTRACT 

We report on the spectroscopic transit of the massive hot-Jupiter CoRoT-Exo-2b observed with the high-precision spectrographs SOPHIE and 
HARPS. By modeling the radial velocity anomaly occurring during the transit due to the Rossiter-McLaughlin (RM) effect, we determine the 
sky-projected angle between the stellar spin and the planetary orbital axis to be close to zero A = 7.2 ± 4.5 deg, and we secure the planetary nature 
of CoRoT-Exo-2b. We discuss the influence of the stellar activity on the RM modeling. Spectral analysis of the parent star from HARPS spectra 
are presented. 
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1. Introduction 

Measurement of the spectroscopic signal during the transit of 
an exoplanet in front of its host star - known as the Rossiter- 
McLaughlin (RM) effect - provides an assessment the trajectory 
of the planet across the stellar disk and, more precisely, the sky- 
projected angle between the planetary orbital axis and the stellar 
rotation axis. This misalignment angle, denoted by A, is a funda- 
mental property of planetary systems that provides clues about 
the process of planet migration. Among the 30 transiting exo- 
planets known so far, A has been reported for only 5 exoplanets 
(HD209458b, Queloz etal.EOOOl HP 18973 3b, Wi nn et al. 120051 
HAT-P-2, Winn et al. 120071 Loeillet et al. 120081 HD149026b, 
Wolf et al. 120071 and TrES-1, Narita et al . 120071 . For all of these 
cases, A is close to zero, as in the solar system, and the stellar 
rotation is prograde relative to the planet orbit. Such measure- 
ments should be extended to other transiting systems to under- 
stand whether this degree of alignment is typical. 

The massive hot-Jupiter CoRoT-Exo-2b (Alonso et al. 
2008) was revealed as planetary candidate by the CoRoT space 
mission (Baglin et al. |2003t and its planetary nature and mass 
was established thanks to ground-based facilities, including 
high-precision spectrographs SOPHIE (Bouchy et al. 2006J and 
HARPS (Mayor et al. |2U03l ). This second CoRoT exoplanet 
is a 3.3 Jupiter-mass planet orbiting an active G7 dwarf star 
(mv=12.6) every 1.743 days. We report here the measurements 
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* Observations made with SOPHIE spectrograph at Observatoire de 
Haute Provence, France (PNP.07A.MOUT) and HARPS spectrograph 
at ESO La Silla Observatory (079.C-0127(F)). The CoRoT space mis- 
sion, launched on December 27th 2006, has been developed and is oper- 
ated by CNES, with the contribution of Austria, Belgium, Brasil, ESA, 
Germany, and Spain. 



of the spectroscopic transit observed with both SOPHIE and 
HARPS spectrographs. These observations were made simul- 
taneously with the space-based photometry with CoRoT. Such 
simultaneous monitoring is useful to assess anomalies in the 
transit parameters due to star spots or transient events. 

Our data permits us to determine the sky-projected angle 
between the stellar spin and the planetary orbital axis, and it 
provids additional constraints on the orbital and physical param- 
eters of the system. Furthermore, our data confirms and secures 
the planetary nature of the transiting body, excluding blending 
of an eclipsing binary with a third star as the cause of the ob- 
served shallow transits. We used HARPS spectra to perform the 
spectroscopic analysis of the parent star. 

2. Observations 

We performed high-precision radial velocity observations of 
CoRoT-Exo-2 (mv=12.6) with the SOPHIE spectrograph, based 
on the 1.93-m OHP telescope (France), and the HARPS spec- 
trograph, based on the 3.6-m ESO telescope (Chile). These two 
instruments are cross-dispersed, fiber-fed, echelle spectrographs 
dedicated to high-precision Doppler measurements based on 
the radial velocity techniques of simultaneous-Thorium calibra- 
tion. SOPHIE was used with its high efficiency mode (spectral 
resolution R=40,000). We reduced HARPS and SOPHIE data 
with the same pipeline based on the cross-correlation techniques 
(Baranne et al. [T9961 Pepe et al. 120021 . We observed CoRoT- 
Exo-2 with SOPHIE on 16 July 2007 and with HARPS on 1 
September 2007. The exposure times were respectively 10 and 
20 minutes on HARPS and SOPHIE corresponding to S/N per 
pixel at 550 nm of 16 and 25, respectively. We obtained the ra- 
dial velocities by weighted cross-correlation with a numerical 
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Table 1. Radial velocity measurements of CoRoT-Exo-2 ob- 
tained by HARPS and SOPHIE during the transit. BJD is the 
Barycentric Julian Date. 



BJD 


RV 


Uncertainty 


-2400000 


[kms- 1 ] 


[kms- 1 ] 


SOPHIE 2007-07-16 


54298.4641 


23.341 


0.026 


54298.4862 


23.285 


0.027 


54298.5030 


23.369 


0.028 


54298.5198 


23.378 


0.027 


54298.5381 


23.123 


0.027 


54298.5550 


22.926 


0.028 


54298.5714 


22.891 


0.029 


54298.5879 


23.023 


0.030 


HARPS 2007-09-01 


54345.5225 


23.371 


0.020 


54345.5298 


23.371 


0.019 


54345.5371 


23.392 


0.018 


54345.5444 


23.360 


0.018 


54345.5517 


23.347 


0.019 


54345.5590 


23.370 


0.018 


54345.5663 


23.456 


0.018 


54345.5736 


23.488 


0.020 


54345.5809 


23.496 


0.022 


54345.5883 


23.394 


0.019 


54345.5956 


23.271 


0.018 


54345.6029 


23.190 


0.018 


54345.6124 


23.060 


0.018 


54345.6197 


23.005 


0.017 


54345.6270 


22.939 


0.018 


54345.6343 


22.928 


0.018 


54345.6417 


23.034 


0.019 


54345.6490 


23.093 


0.019 


54345.6563 


23.107 


0.019 


54345.6636 


23.090 


0.019 


54345.6709 


23.034 


0.020 



G2 mask constructed from the Sun spectrum atlas including up 
to 3645 lines. We eliminated the first 8 blue spectral orders con- 
taining only noise. Radial velocities are given in Table Q] and 
displayed in Fig.Q] 

3. Rossiter-McLaughlin modeling 

The RM effect corresponds to a distortion of the spectral lines 
observed during a planetary transit due to stellar rotation. The 
transiting body hides some of the velocity components that usu- 
ally contribute to line broadening resulting in an Doppler-shift 
anomaly (see Otha et al. 2005 , Gimenez et al. |2006b ; Gaudi & 
Winn |2007] l. 

To model this RM effect, we used the analytical approach 
developed by Otha et al. (1200 5). The complete model has 12 
parameters: the orbital period P; the mid-transit time T c ; the ec- 
centricity e; the angle between the node and periastron co; the 
RV semi-amplitude K; the velocity zero point Vb (these first six 
are the standard orbital parameters); the radius ratio r p /R s ; the 
orbital semi-major axis to stellar radius a/R s (constrained by 
the transit duration); the sky-projected angle between the stel- 
lar spin axis and the planetary orbital axis A; the sky-projected 
stellar rotational velocity v sin /; the orbital inclination z; and the 
stellar limb-darkening coefficient e. For our purpose, we started 
with the orbital parameters and photometric transit parameters 
as derived by Alonso et al. (2008). We fixed the linear limb- 
darkening coefficient e=0.78, based on Claret (2004) tables for 
filter g' and for the stellar parameters derived in Sect. [4] Our free 



parameters are then A and v sin /. We introduced two additional 
parameters: the offset velocity of HARPS and SOPHIE, Aharps 
and Asophie, which differ from Vq due to the stellar activity. 
We determined the vsinl independently from SOPHIE cross- 
correlation functions (CCFs) to be 9.5 ±1.0 kms -1 and from 
HARPS CCFs to be 10.7 + 0.5 kms 1 with the calibration tech- 
niques described by Santos et al. (2002). However, we decided 
to leave it as free parameter in our fit. 

The result of our fit, displayed in Fig.[T]and listed in Table|2] 
first shows that the stellar rotation is prograde relative to the 
planet orbit. During the first part of the transit the starlight is 
redshifted, indicating that the planet is in front of the approach- 
ing (blueshifted) half of the stellar disk. During the second part 
of transit, the sign is reversed as the planet moves to the receding 
(redshifted) half of the stellar disk. The sky-projected angle be- 
tween the stellar spin axis and the planetary orbital axis A is close 
to zero. The projected rotation velocity of the star v sin / deter- 
mined by our RM fit (11.85+0.5 kms -1 ) seems slightly larger 
than our spectroscopic determination (2-cr greater). Previous 
studies by Winn et al. ([2005) showed that the v sin / measured 
with Otha formulae was biased toward larger values by ap- 
proximatively 10%. But, as already suggested by Loeillet et al. 
(2008), it may be due to the differential rotation of the star from 
equator to pole. Considering the exoplanet crosses the star near 
its equatorial plan, the fitted v sin / corresponds to the maximum 
value. Note that if we fix v sin / at the spectroscopic value, it does 
not change the value of the fitted A angle. 

We made the 2 epochs of RM observations at a minimum 
stellar flux (see Fig. 1 of Alonso et al. 120081) . indicating that 
the stellar spots were at their maximum phase of visibility. 
Following the Saar & Donahue (1997) relation giving the ex- 
pected RV jitter as a function of v sin / and spot filling factor, we 
found that CoRoT-Exo-2 is expected to present RV variations of 
up to 200 ms -1 peak-to-peak with a period of 4. 5 days. The stan- 
dard deviation of RV residuals (56 ms ) found by Alonso et al. 
(2008) is in agreement with this value. Such an activity-related 
RV variation should then change locally the apparent slope in the 
RV orbital curve. The maximum effect occurs at the maximum 
phase of stellar spot visibility, and should induce an apparent in- 
crease in the semi-amplitude K of up to 40 ms" . This explains 
why our fit in Fig.[TJis not perfect outside of the transit. If we in- 
crease K in our fit or leave it as a free parameter, it significantly 
improves the fit and slightly decreases the value of v sin / and A 
(see Table|2]i. 

We also did a combined fitting of the photometry and the 
whole set of RV measurements. On each of the out-of-transit 
measurements, we inserted an additional error on the RV data to 
take the stellar activity into account. We chose this value as 56 
ms -1 , corresponding to the standard deviation found by Alonso 
et al ( 120081) . This correction is justified since the action of ac- 
tivity on the points taken at random out-of-transit phases can be 
assumed as random for these points, while during transit we have 
sets of points with the same activity level throughout. The fitting 
was done using a Markov Chain Monte Carlo (MCMC) with 
a Metropolis-Hastings Algorithm for the decision process. We 
used the models of Gimenez 2006a and 2006b for photometry 
and spectroscopic transits respectively. A quadratic law of limb 
darkening was used. For the photometry, we used the fitted pa- 
rameters found by Alonso et al. (2008). For the spectroscopy pa- 
rameters, we chose the V-band, from tables published by Claret 
2000 for the stellar parameters derived in Sect. [4](m + = 0.748, 
m_ = 0.256). The MCMC was performed over 20000 accepted 
steps after 5 000 steps of a burn-in period. The result of the com- 
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Fig. 1. Phase-folded radial velocity measurements of CoRoT-Exo-2 during the transit of the planet with SOPHIE (dark circle) 
and HARPS (open circle). The solid line corresponds to the Rossiter-McLaughlin model ajusted to these data assuming the semi- 
amplitude K-563 ms _1 from Alonso et al. (2008). The dotted line corresponds to the Rossiter-McLaughlin model with K as free 
parameters. 



Table 2. System parameters of CoRoT-Exo-2. The reduced y 1 
was computed assuming 24 degrees of freedom. 

Fixed parameters from Alonso et al. 2008 



p 


1.7429964 days 


Tc 


54237.53562 


e 


0.0 


V 


23.245 kms 1 


r p IR s 


0.1667 


a/R s 


6.70 


i 


87.84 deg 


e 


0.78 (from Claret) 


Adjusted parameters with X"=563 m s 1 


v sin / 


11.85± 0.50 kms" 1 


A 


7.2+4.5 deg 


Sharps 


-21.5 ± 5 ms-' 


As OPHIE 


+21.5 ± 12 ms" 1 


reduced x 1 


1.43 


Adjusted parameters with K as free parameter 


K 


656±27 ms H 


v sin / 


11.25+ 0.45 kms" 1 


A 


5.0+4.0 deg 


Sharps 


-25.0 + 4.5 ms" 1 


&s OPHIE 


+25.5 ± 11 ms' 1 


reduced x 1 


1.01 



Combined MCMC fit 



K 

vsin/ 
A 

Sharps 

As OPHIE 

reduced \ 2 



613+14 ms" 1 
11.46+°-^ kms" 1 
7.1+5.0 deg 
-22.5 + 4.5 ms" 1 
+23.5 ± 11 ms" 1 
1.10 



bined fit is presented in Tabled and is in full agreement with the 
other approachs. 

The cross-correlation function (CCF) corresponds more or 
less to an average of all the spectral lines (see top of Fig. In 



order to characterize the behavior of the spectral lines during the 
transit, we computed the difference between the HARPS CCFs 
corrected from the orbital velocity and a reference CCF taken out 
of the transit (more exactly an average of the 3 first exposures). 
This difference was computed at 5 epochs identified and labeled 
in Fig. Q] : (1) just before the ingress, (2) maximum of the RM 
effect, (3) mid-transit epoch, (4) minimum of the RM effect, (5) 
just after egress. These differences ACCF = CCFref-CCF# are 
displayed in Fig. [2] and clearly show the spectroscopic anomaly 
shifting from the blue side (2) to the red side (4) of the CCF. 
During the transit, the depth or contrast of the CCFs is system- 
atically larger, reflecting the renormalization effect of the CCF, 
which maintains a constant surface. 



The observation of the spectroscopic transit of CoRoT-Exo- 
2b allows us to confirm definitively that the transiting candidate 
provided by CoRoT occurred at the central star (and not at a 
background star inside the CoRoT PSF). Furthermore, if we as- 
sume that the system is not diluted by an other star inside the 
HARPS or SOPHIE PSF, the RM anomaly reveals that the tran- 
siting body has a planetary size (from the RM anomaly ampli- 
tude) and planetary mass (from the RV slope outside the tran- 
sit). In the case of an eclipsing binary whose light is diluted 
with a brighter third star, one should assume that the spectral 
lines of the fainter eclipsing binary move relative to the lines 
of the bright star and thus change the blended line-profiles. In 
such a configuration, one should consider not only the flux ratio 
but the v sin /, velocity zero point, and spectral type of the two 
systems. In our present case, we did not find a configuration of 
a blended eclipsing binary that could simultaneously reproduce 
the RV anomaly and the photometric light curve. Furthermore, 
we computed RVs using different cross-correlation mask with- 
out significant changes in the shape and amplitude of the RM 
anomaly. We, thus, consider that the spectroscopic transit con- 
firms and secures the planetary nature of the transiting body. 
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Fig. 2. (Top) Averaged cross-correlation function of CoRoT- 
Exo-2. (Bottom) Cross-correlation differences computed at 5 
different epochs (see text) illustrating the behavior of the spectral 
lines during the transit. 



4. Spectroscopic analysis of CoRoT-Exo-2 

We performed the spectroscopic analysis of the parent star using 
the HARPS spectra. We corrected individual spectra from the 
stellar velocity, rebinned to a constant wavelength step of 0.02 
A, and co-added spectral order per spectral order giving a S/N 
per pixel at 550 nm of about 80. We determined the effective 
temperature first from the analysis of the Ha line wings, provid- 
ing a temperature of 5450+120 K. In spite of the quite low S/N 
of the combined spectra, it appears that the star is at the border of 
the temperature domain in which the Ha line wings are a good 
temperature indicator (from 5500 to 8500 K). We then checked 
this result with other methods. We performed synthetic spectra 
fitting using LTE MARCS atmosphere models (Gustafsson et al. 
2005), which are well adapted for this range of temperature. We 
compared the synthetic spectra, previously convolved by the in- 
strumental profile and a rotational profile with the v sin / value 
previously measured, to the observed one. The best-fit model 
yields a slightly higher temperature, but is still in agreement with 
the Ha estimate. Another analysis, using equivalent width mea- 
sures of Fel and Fell lines, was also carried out and yields simi- 
lar results. The adopted stellar parameters are T e fr =5625+ 120 K, 
logg=4.3+0.2 and [M/H]=0.0±0.1, which correspond to a G7V 
type star with a solar metallicity. With these values, we derived 
the star's luminosity and mass with StarEvolv stellar evolution 
models (Siess 2006; Palacios, private communication). We com- 

1 /3 

bined these estimates of the star's mass to the M s /R s value 
provided by the light curve analysis to derive the final star's 
mass and radius values in a consistent way between spectro- 
scopic and photometric analyses. The method allows us to get 
rid of the large uncertainty that affects the estimate of the gravity 



and to take advantage of the excellent quality of the light curve. 
The method will be detailed in a forthcoming paper devoted to 
the fundamental parameters of the first CoRoT planet host stars, 
based on UVES spectra. The adopted stellar mass is 0.97+0.06 
M Q and the stellar radius is 0.90+0.02 R . Interestingly, the 
solar-like metallicity of the parent star and large radius of the 
planet is consistent with the trend that heavy element content 
in the planet and stellar metallicity are correlated (Guillot et al. 
2006). According to stellar evolution models (Lebreton, private 
communication), the age of the star could be between 0.2 and 4 
Gyr if the star is on the main sequence. However, the presence of 
the Li I absorption line and the strong emission line core in the 
Call H and K lines, suggest that the star is still close the ZAMS 
and could be thus younger than 0.5 Gyr in full agreement with 
the observed stellar activity and the measured rotation period. 

The knowledge of the main rotational period of CoRoT-Exo- 
2 determined from the light curve (4.54 days) and the spectro- 
scopic v sin / determined from HARPS and SOPHIE CCFs (10.3 
kms" 1 ) may be used to independently estimate the minimum 
radius of R s sin/= 0.92 R Q in very good agreement with our 
previous determination based on spectral classification. We note 
that this estimate, based on the well-determined stellar rotation 
thanks to the high-precision CoRoT light curve, does not depend 
on any spectral classification. On the other hand, if we assume 
the stellar radius from spectral analysis, we can deduce that sin / 
is close to 1, indicating a further constraint on the alignment of 
orbit and stellar spin. 

5. Conclusions 

In addition to the previous 5 transiting exoplanets where A angle 
have been reported, CoRoT-Exo-2b presents a prograde orbit rel- 
ative to the stellar rotation and an angle A = 7.2 + 4.5 deg, close 
to zero. Our observations illustrate and demonstrate the capabil- 
ity of extending the reach of the RM technique to relatively-faint 
host stars (mv>12) like the CoRoT targets even with a 2-m class 
telescope. 
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